Direct Evidence for a Two-component Tunnelling Mechanism 
in the Multicomponent Glasses at Low Temperatures 
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The dielectric anomalies of window-type glasses at low temperatures (T < 1 K) are rather suc- 
cessfully explained by the two-level systems (2LS) tunneling model (TM). However, the magnetic 
effects discovered in the multisilicate glasses in recent times an d also some older data from 

mixed (Si02)i-;c (K20)a; and (Si02)i-a;(Na20)a; glasses indicate the need for a suitable general- 
ization of the 2LS TM. We show that, not only for the magnetic effects 0,0] but also for the mixed 
glasses in the absence of a field, the right extension of the 2LS TM is provided by the (anomalous) 
multilevel tunneling systems approach proposed by one of us. It appears that new 2LS develop via 
dilution near the hull of the SiCU-percolating clusters in the mixed glasses. 

PACS numbers: 61.43.Fs, 07.20.Mc, 61.43.Hv, 77.22.-d, 77.84.Lf, 65.60.+a 
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Tunneling systems are useful for many theoretical in- 
vestigations in condensed-matter physics: from laser- 
matter interaction, to amorphous solids and disordered 
crystals, to hydrogen storage in metals. For their sim- 
plicity, only the two-level systems (2LS) are usually em- 
ployed instead of the more realistic multilevel systems; in 
this Letter a situation will be discussed where the use of 
a number of states greater than two is essential. More- 
over, new insight will be given on the role of percolation 
and fractal theory in the TM of multicomponent glasses. 

Glasses at low temperatures are interesting for their 
rather universal physical properties attributed to the low- 
energy excitations characterising all kinds of amorphous 
solids. The 2LS tunneling model (TM) has been suc- 
cessful in the semi-quantitative explanation of a variety 
of interesting thermal, dielectric and acoustic anomalies 
of structural glasses at temperatures T < IK [6]. The 
physics of cold glasses is important also because of its 
link with the mechanism of the glass freezing transition 
0- Beside the linear in T behavior of the heat capacity 
C p , it is believed that the linear in ± lnT behavior of the 
real-part of the frequency-dependent dielectric constant 
e'(u>) represents a cogent characterization of the glassy 
state at low temperatures. 

Fig. 1 (inset) shows the behavior of the T-dependent 
part of e', Ae'/e' = [e'(T) - e'(T )]/e'(T ), (where T (uj) 
is a characteristic minimum) for vitreous SiC>2. It can 
be seen that linear regimes in — In T for T < To and 
+ lnT for T > Tq are observed, and roughly with slopes 
S- = —2S and S+ = S > 0, or in a -2:1 ratio. According 
to the 2LS TM, in fact, we have the expressions 
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where we neglect, for low u>, the frequency-dependence in 
the RES part, z m i n ,max = &0min,max/2kBT and where 
r is the phenomenological 2LS relaxation time given by, 
with E = 2k B Tz 
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FIG. 1: Dielectric signature of pure a-SiC>2 (inset) and AlBa- 
SiO (main) glasses. Si02 data [g], fitted with Eq. {T}, display 
a -2:1 2LS TM behavior. AlBaSiO data [J display rather a 
-1:1 behavior, yet could be fitted with Eq. (JTJ (dashed line) 
[j| with a large A um i„ =12.2 mK 2LS tunneling parameter. 
We have fitted all data with a more realistic Aomin =3.9 mK 
and best fit parameters from Table 1 using Eq.s |T} and ([5]) 
(driving frequency u> =1 kHz). 

Here, Ao is the tunneling parameter of a 2LS hav- 
ing energy gap E, A 0m in and A Qmax are its phe- 
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nomenological bounds, 7 is an elastic material param- 



eter of the solid and r. 
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P is the probability per unit 



volume and energy that a 2LS occurs in the solid and -p\ 
is the average square 2LS electric dipole moment. 

Indeed, from expressions {J) we see: 1) The so-called 
resonant (RES) contribution has the leading behavior 
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2) the relaxational (REL) contribution has, instead, the 
leading behavior 
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Thus, the sum of the two contributions has a V- 
shaped form, in a semilogarithmic plot, with the mini- 
mum occurring at a T roughly given by the condition 
VT min (k B T) ~ 1, or k B T (uj) ~ (^w) 1 / 3 . e e r is here 
the bulk of the solid's dielectric constant and we see that 
a -2:1 characteristic behavior is justified by the TM with 
the T > T slope given by S = P-p\jZt^r- 

This behavior is observed in pure a-Si02 [H (with the 
parameters of Table 1, x = 0, from our own best fit to 
Eq. |T])), however in most multicomponent glasses it is 
rather a V-shaped curve with a (roughly) -1:1 slope ratio 
that is often observed. Fig. 1 (main) shows this phe- 
nomenon for the multisilicate glass A^Os-BaO-SiC^ (in 
short, AlBaSiO), which has been extensively investigated 
in recent times due to its unexpected magnetic field re- 
sponse [J Also, Fig. 2 shows the behavior of the 
dielectric constant vs T for the glasses of composition 
(Si02)i-x(K20) 2; containing a molar concentration x of 
alkali oxide |4j. It is seen that a S-/S+ slope ratio of 
roughly -1:1 is observed, with the slope definitely chang- 
ing with x (and faster for T > Tq). These data, thus far 
unexplained by the 2LS TM, call for an extension of the 
accepted TM and we show in this Letter that a simple 
explanation can be given in terms of the very same new 
(anomalous) tunneling systems (ATS) that have been ad- 
vocated by one of us to explain the magnetic response of 
AlBaSiO and other multicomponent glasses d, In 
view of the interest for these materials in low-T ther- 
mometry, and on fundamental grounds, such explanation 
appears overdue to us. Moreover, "additional" TS (be- 
side the standard 2LS) of the type here advocated were 
already called for in [J] and earlier papers. 

In a rather general fashion, the TS can be thought 
of as arising from the shape of the theoretical energy- 
landscape E({ri}) of a glass as T is lowered below the 
glass freezing transition Tf. Many local and global min- 
ima develop in E({ri}) as T — > 0, the lowest-energy min- 
ima of interest being made up of n w = 2,3,... local wells 



separated by shallow barriers. These local multiwelled 
potentials are our TS and it seems reasonable that the 
n w = 2 - welled potentials will be ubiquitous in this pic- 
ture. These should be thought of as effective represen- 
tations of local "tremblements" of the equilibrium posi- 
tions {r) '} of some of the glass ions' positions (unlike 
in the disordered crystal case, where the TS ought to be 
rather well- localized dynamical entities). Hence, just as 
the n w — 2 - welled case is possible, so ought to be the 
n w =3,4,...- welled situations which would also be lo- 
cal rearrangements involving a few atoms/ions. So long 
as their energy parameters obey the usual uniform distri- 
bution advocated by the TM, however, most of these n w - 
welled potentials should present the very same physics as 
the n w — 2 cases and thus in practice the n w distribution 
cannot be resolved experimentally in a pure glass. 
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FIG. 2: Dielectric signature of mixed (SiO^i-^K^O)^ 
glasses as function of T and x [J]. Fitting parameters from 
Table 1 using Eq. ([l]) and ([5]) from our theory (driving fre- 
quency ui —10 kHz). 

All change if the glass is made up by a mixture 
of network-forming (NF) ions (like those of the SiC>4 
or (AIO4) - tetrahedral groups) as well as of network- 
modifying (NM) ions (like K+ or Na+, or Ba 2+ , Pb 2+ , 
...) which, these last ones, could act as nucleating cen- 
tres for a partial devitrification of the glass, as is know 
to occur in the multicomponent materials 1^, U | ■ Simu- 
lations and experiments in the multisilicates have shown 
that the NM-species in part destroy the networking ca- 
pacity of the NF-ions and form their own pockets and 



channels inside the NF-network [12|. Hence, n w > 2 mul- 
tiwelled systems inside these NM-pockets and -channels 
should follow some new energy-parameters' distribution 
forms when some degree of devitrification occurs, leading 
to entirely new physics. 

One of the present Authors has proposed that precisely 
this situation occurs inside the magnetic-sensitive multi- 
component glasses [3, 5] , and in this Letter we show how 
this theory explains the B = dielectric data of Fig.s 
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1-3 as well. Using the theory of 0, @| we arrive at the 
expression, for the contribution to the dielectric anomaly 
from the advocated ATS: 
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Here we have again neglected, for low-w, the frequency- 
dependence in the RES part, we put y = D/D m i n and 
TAmax is the largest phenomenological ATS relaxation 
time given by 
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Moreover D is the lowest-lying energy gap of a multi- 
level ATS, r is an appropriate elastic constant and P* 
is the (slightly renormalised) probability per unit vol- 
ume that an_ATS occurs in the NM pockets and chan- 
nels, with p\ the average square ATS dipole moment. 
This description is linked to a distribution function 0, 0] 
for these ATS favoring near-degenerate energy gaps D 
bound from below by D m i n . In turn, this produces 
an overall density of states g(E) = g 2 LS + 9ats{E) ~ 
2P + (2nP*/E)9(E - D min ) that is roughly of the form 
advocated in |4j and by some other earlier Authors (e.g. 
[l3T |) to explain anomalies not accounted for by the stan- 
dard 2LS TM. 

Manipulation of the expressions in shows that: 1) 
The RES contribution from the ATS has the leading be- 
havior (for T < D m i n /2k B , e'\ARES is a constant) 
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2) the REL contribution is, instead, characterised by the 
leading form 
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Thus, the V-shaped semilogarithmic curve is somewhat 
lost. However adding the 2LS (Eq. {TJ) and ATS (Eq. 
([5])) contributions together one does recover a V-shaped 
curve with a slope S~ ~ — 2S basically unchanged for 
T < Tq and an augmented slope S+ = S + Sats f° r 
T > T with Sats = 7TrP*pJ/6eoe r k B T that for T < 
Dmin/kB may approach 2S and thus (qualitatively) a 
-1:1 slope ratio. 
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FIG. 3: Dielectric signature of mixed (Si02)i-a;(K20)a: 
glasses as function of T and ui for x =0.2 Fitting pa- 
rameters from Table 1 using Eq. |T} and ((5JI from our theory. 



glass 
type 


X A2LS 7 AaTS Dmin T 

1CT 5 10" 8 sJ 3 10" 5 K 10" 6 sK 5 


SiOa 


47.2 5.30 


AlBaSiO 


- 116.2 13.40 264.7 0.65 69.73 


K-Si 


0.05 104.1 1.33 75.5 0.87 3.55 


K-Si 


0.08 146.5 1.23 130.0 0.87 3.97 


K-Si 


0.10 158.5 1.15 160.0 0.87 5.08 


K-Si 


0.20 239.5 0.82 281.9 0.87 6.44 



TABLE I: Extracted parameters for the glasses; K-Si stands 
for the (Si02)i-x(K20)a; glasses. In all of the best fits we 
have employed the values A Um in =3.9 mK and Ao, nal =10 K 
extracted from fitting the pure SiC>2 data of Fig. 1 (inset). 



We have fitted expressions |T]) and ([5]) to the data for 
AlBaSiO in Fig. 1 (main) and to the s-dependent data 
for (Si02)i- 2; (K20) I in Fig.s 2 and 3, obtaining in all 
cases very good agreement between theory and experi- 
ments [14] . Fig. 3 shows the fit of our theory to the 
frequency-dependent data for x = 0.2. In all these best 
fits we have kept the value of Ao m i„ = 3.9 mK fixed, as 
obtained from our pure Si02 fit, and the value of D m i n 
also independent of x and w. The idea is that these 
parameters are rather local ones and should not be influ- 
enced by NF/NM dilution. Table 1 reports all the (2LS 
and ATS) parameters used for our best fits and for Al- 
BaSiO the parameters differ somewhat from the values 
A ATS = 1.42xl(T 3 and D min = 0.03 K used in 5]. This 
shows that a more microscopic version of this model is 
called for. Fig. 4 shows the dependence of the prefactors 
with x. It can be seen that, as expected, the ATS prefac- 
tor A ats = kP*Px/ eQ£rD m i n scales linearly with x, an 
excellent confirmation that the "additional" TS of [1, [l3| 
are indeed our ATS forming near the microcrystallites 
within the NM-pockets and channels. It can be seen, in- 



4 



stead, that the 2LS prefactor A^ls = PPo/ e o^r of our 
fits also increases, though less rapidly, with x (a decrease 
like 1 — x would be expected). We propose (adopting 
a NF/NM percolation picture) that new, "induced" 2LS 
form with alkali dilution near the NF/NM surface of the 
NF percolating clusters as x is increased from 0. This 
leads to the expression Ab u ik{l — x)+A sur fP(x)x^ for the 
2LS prefactor, with A(, u ik, A sur f and / fitting parameters 
and P{x) the percolation probability function (P(x) ~ 1 
for small x). Our best fit leads to / = 0.81, in rather good 
agreement with the euristic expression / = 1 — (D — D s )v 
(D is the fractal dimension of the percolating cluster, 
D s < D of its "elastic" surface (not necessarily the hull) , 
v the connectedness length's exponent) one would infer 
from elementary fractal or percolation theory. 




A 2LS = A B (1-X)+A S P(X)X 
A„=47.15, A„=748.07, f=0.81 
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FIG. 4: The 2LS and ATS prefactor parameters (xlO 5 ) for 
all glasses (from Table 1) as a function of x. Our data fit well 
with our theoretical expectations (full lines). 

Our theory, with parameters from Table 1, can also ac- 
count for the T— and x— dependence of the heat capacity 
C p (T,x) (data from [4]; for AlBaSiO see also @). 

In summary, we have shown that there is direct evi- 
dence in zero magnetic field already for the multiwcllcd 
ATS advocated to explain the magnetic field effect in the 
multicomponent glasses. The relevance of multiwelled 
tunneling systems in the multicomponent glasses is a new 
and unexpected finding in this field of research. Our 
work predicts that the magnetic response of the alkali 
silicate glasses should be important and scale like the 
molar alkali concentration x. At the same time the -1:1 
slope ratio problem of the standard TM has been given a 
simple explanation in terms of our two-component tun- 
neling model. For AlBaSiO, using the value of P* ex- 
tracted from C p [3] and that of Aats given in Table 1, 
we estimate the value of the electic dipole of the effective 

anomalous tunneling entity: p e // = \fp\ ~0.4 D. This 
small value confirms 0, Q that this physics must come 
from the coherent tunneling of a small ionic cluster (the 
very same origin of the large value of D m i n ). Details of 



our analysis will be published elsewhere. 

The present work suggests how new experiments 
on mixed glasses could help in deciding whether our 
two-component tunneling model 0, Q or the nu- 
clear electric quadrupole-moment approach [l5[ is the 
correct explanation for the magnetic-field effects in 
the cold glasses. Amorphous solids of composition 
(Si02)i- a; (MO) 2 ; should be investigated for different con- 
centrations x of NM-species at ultralow temperatures and 
in the presence of a magnetic field. If M is a divalent 
metal (ideally Ca, but also Mg, Sr and Ba) with no or 
a tiny percentage of stable isotopes possessing a nonzero 
nuclear quadrupole moment, then a magnetic response 
scaling like x would single out the present approach. The 
system (Si02)i- :E (CaO) ;r (with a small percent of PbO 
favoring nanocrystal formation) appears to be the ideal 
candidate. We remark, incidentally, that in Q the amor- 
phous system Si02+ 2 :Cj / H z , devoid of nuclear quadrupole 
moments, displayed magnetic effects quantitatively sim- 
ilar to those found in glasses with nuclei carrying such 
moments [l| at ultralow temperatures. 
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